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Dispersion of "Dispersionless Band": 
Comments on L. Brey and H.A. Fertig paper Electronic States of Graphene 
Nanoribbons Studied With The Dirac Equation 

Lyuba Malysheva and Alexander OnipkcQ 
Bogolyubov Institute for Theoretical Physics, 03680 Kyiv, Ukraine 
(Dated: May 8, 2008) 

As a particular application of the earlier proposed model of graphene as a macromolecule, we 
have found the exact analytical expression of dispersion relation for the band of edge states in 
graphene zigzag ribbons. This band is often referred to as "dispersionless band" or "zero modes". 
The obtained results contrast conclusions regarding edge states given in the referenced paper and 
show where the earlier given description is valid and where it is not correct. 



INTRODUCTION 

Dramatic effects of zigzag edges on the band struc- 
ture and density of states of finite-size layers of 2D 
graphite were demonstrated by first-principles calcula- 
tions in 1993 Tj. In many subsequent studies, pe- 
culiar graphene properties associated with edge states, 
have been studied in many details by different methods 
0, 0, 0, H, @, 0, H S E3] • In particular, exponential behav- 
ior of the density of states and dispersion near the Fermi 
energy has attracted much of attention. An exponential 
decrease of edge state energy observed in numerical stud- 
ies [l|, Q , was explained analytically by Brey and Fertig 
Q, henceforth referenced as BF. In this paper, electron 
states of graphene ribbons (GRs) were studied with the 
use of the Dirac equation, that is within a rclativistic 
model [llj , originating from the tight binding descrip- 
tion of infinite 2D graphite [l2|, [13]. To take into 
account edge effects, certain boundary conditions have 
been imposed on the wave function, separately, for A and 
B sublattices indicated in Fig. la. An impressive quanti- 
tative agreement between tight-binding calculations and 
solutions of the Dirac equations was demonstrated by 
a number of illustrative examples. However, from the 
methodological point of view, such an approach suffers 
from certain weaknesses, that will be addressed here later 
on. 

As mentioned, the edge-state spectrum of zigzag GRs 
was quantified by BF in an analytical form. It relates 
the decrease of energy of edge states with an increase of 
separation of longitudinal wave vector ky F from a zero- 
energy point of 2D graphite dispersion relation. One of 
these so-called K points is shown in Fig. 1. According to 
BF, the spectrum of edge states is described by two ap- 
proximate equations, reproduced below by Eqs. ([I]) and 
<[2|) in slightly different notations. 

Here we show that this result is correct only partly. 
In fact, for the larger part of edge states, the edge-state 
energy is governed by an exponential behavior with dif- 
ferent pace and crossover from one exponential depen- 
dence for smaller wave vectors to another for larger wave 
vectors. These dependencies, represented in Eqs. Q and 



(fTTT) . follow from the exact treatment of graphene as a 
macromolecule [IH . Our equation ((4J agrees but Eq. (fTT|) 
disagrees with the earlier prediction. We see the reason 
of the divergence in that the BF approximate descrip- 
tion is not good enough for zigzag edges and therefore it 
does not catch all peculiarities of electronic properties in 
a very narrow region near the Fermi energy, where the 
edge states occur. 



HOW DOES ELECTRON ENERGY GO TO ZERO 
IN ZIGZAG GRAPHENE RIBBONS? 

To begin with, let us write the well known lin- 
ear dispersion relation for an infinite graphene sheet 

± ^T to \A fc x F ) 2 + i k y F ) 2 : where t is the hopping inte- 
gral between nearest-neighbor carbons in the honeycomb 
lattice. Label BF indicates that the reference point in k- 
space (K point of the 2D graphite band structure) is the 
same as in the referenced paper. According to BF, for 
a graphite sheet of finite width L — \f?>a{N — 1/3) and 
an infinite length (a zigzag GR as it appears in Fig. la) 
the dispersion relation changes. Namely, if longitudinal 
component fc^ F exceeds the critical value ky = l/L, the 
transverse wave-vector component becomes imaginary, 

L.BF _ Ah 



4- A 
±— ta 



fc.!? F > kl 



and satisfies equation 



fccoth(fcL) = k* F . 



(1) 



(2) 



This is a slight modification of Eq. (6) as it appears in 
Ref. (|. 



BF 



> 



It is easy to verify that under the condition ky 
2ky, the dependence of real solution to Eq. ||3J| on ky F is 
accurately (with error less than 1%) reproduced by 



k = k 



BF/ 



. s '(l-2e-". 
The use of this result in Eq. (p} yields 



(3) 



e = ±V3taky F exp{~k^ F L) = ±V3tqexp{-V3Nq), 



(4) 



2 



which finalizes a fully analytical description of the edge- 
state spectrum. Here and henceforth, N is assumed to 
be a large number in all approximate relations. The sec- 
ond equality represents a new notation, ak^ F = q, which 
is more suitable for further consideration. A point to 
note is that the divergence between Eq. ^ and ener- 
gies calculated from Eqs. ([I]) and ([2]) is noticeable only 
for F L < 2, see inset in Fig. 2. Common limitations of 
these two approximations will be clarified via comparison 
with the exact results discussed next. 



nus branch of Eq. ([5]) can be rewritten as 
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FIG. 1: Zigzag graphene ribbon of width L = v3a(A r — 1/3) 
and infinite length, represented by (a) section of honeycomb 
lattice, periodic in A and B sites (2D graphite model), and (b) 
by N C-C coupled polyacene chains (macromolecule model); 
a is the minimal translation distance. Cross on k y axis in- 
dicates zero-energy point of graphene dispersion relation ((5|. 
Counted from this point, wave vector (k x ,q/a) is equivalent 
to (fc^ F , fcy F ). Quarter of the hexagonal Brillouin zone of 2D 
graphite is also shown. 

The edge-state spectrum can be treated with the use of 
the model shown in Fig. lb. In the framework of canoni- 
cal tight-binding description, the band structure of zigzag 
GRs is determined by [15[ 



±t 



x y^l ± 4 cos(ak y /2) cos(VSak x /2) + 4 cos 2 (ak y / 2), 

(5) 

where the upperscript ± refers to the sign plus or minus 
under the root, k y has the meaning of the longitudinal 
wave vector, < ak y < 7r, and values of ak x = k/\/3 
should be found from equation 



sin k N 



=F2 cos(ak y /2). 



(6) 



sinK±(iV + 1/2) t— y^y 
Note that near zero-energy point, k = ^-(0, |), the mi 



E~ = ± 



V3 



t^(q = ak y -2TT/3y 2 + ( K -/3) 2 , (7) 



making obvious the correspondence between BF's and 
our notations: |afc^ F | = \ak x \ = «~/\/3, and ak^ F = 
q = aky — 2ir/3. 

For the minus branch of dispersion relation ([5]) , Eq. (O 
supports only real solutions n~ , v— 0,1,2,. ..,N—1, if 0< 
aky <27r/3+<f , q°= 2arccos[iV/(2iV+l)]-27r/3. However, 
if q>q c , Eq. © has N— 1 real and one imaginary solu- 
tion. The latter is denoted below as Kq . For N>>1, 
the critical value of the wave vector, q c —(\^3N)^ 1 <<l, 
coincides with aky in the BF theory. 

It is the smallest of fc^-dependent solutions k~ that be- 
comes imaginary, n^lky) — > iS(q). In this case, electron 
energies are bound to the interval \E\/t < (2N)~ 1 , where 



E~ = ±t 



sinh(<5/2) 
sinh(iV + 1/2)5' 



(8) 



and S satisfies 

sinh NS 



smh(iV + l/2)<5 



= cos(g/2) - V3sin(q/2). (9) 



Henceforth, this energy interval is referred to as the band 
of edge states. Distinct from Eqs. {J) and ([2]), Eqs. §8§ 
and ((9]) describing this band are exact. The dispersion 
of edge states predicted by the exact and approximate 
equations is represented in Fig. 2. 

For q c < q « 1, the band of edge states is repro- 
duced by Eqs. dJ)-© reasonably well. Equation ((4]) 
gives nearly equally good description excluding the inter- 
val 1 < q/q c < 2. However, within the larger part of the 
actual interval q c < q < 7r/3, the divergence between the 
approximate and exact descriptions is dramatic. Here, 
we mean the difference between functional forms e(q) and 
E~ (q) rather than the divergence in numbers. 

To specify an analytical expression of edge-state en- 
ergy for larger values of q, one can solve Eq. ([9]) in the 
same approximation that leads to Eq. ([!]). For N5 >> 1, 
the solution to this equation can be expressed as 5 = 
—2 In [cos(<7/2) — \/3 sin(g/2)] . In fact, this approximate 
expression gives fairly accurate results for N5 > 2 and 
can be rewritten in a more compact form 



-2 In [2sm[(7r/3-g)/2]] : 



(10) 



showing that when the wave vector approaches its max- 
imal value, q — > 7r/3, the electron (hole) energy goes to 
zero as 



±*(tt/3 - q) 



2N 



(11) 



In contrast, according to Eqs. I©-©, the edge-state en- 
ergy has though small but finite energy at q = 7r/3. In 
general, the obtained dependence S(q) = V3k(ak® F =q) 



3 



is pronouncedly different from its analogue in the BF de- 
scription. Nevertheless, Eq. ([4]) can be retrieved from 
Eqs. © and (TO]) under the condition q c « q « 1, 
where the restriction on q from above sets the limit of 
applicability of Eq. (6) in Ref. [|| . 

Despite rather crude approximations made under the 
passage from Eqs. (jHJ) and ([9]) to Eq. (fTTj) . the latter works 
reasonably well even for q as 0.3, see Fig. 2. The com- 
parison of the exact, Eq. |5]), and approximate, Eq. (g]) 
and Eq. (1 1 1 [) . dependencies in this figure illustrates two 
regions of the wave vector, where electron/hole disper- 
sion is qualitatively different, although it is exponential 
in both regions. The crossover from one exponential be- 
havior to another occurs within an interval that can be 
approximately designated as it/ 12 < q < tt/6. This spe- 
cific feature of edge-state spectrum is totally lost in the 
BF description, illustrated by the curve (4) in Fig 2. 



E It 




FIG. 2: Energy of electron edge states in a N = 10 wide zigzag 
nanoribbon (semi-logarithmic scale). Exact and approximate 
dependencies on the longitudinal wave vector are marked by 
corresponding equations. Curves calculated from (JTJ, @ and, 
alternatively, from ((4]) are undistinguishable in this scale. In- 
set: same data for small values of q in non-logarithmic scale. 

In connection with the above discussion, frequently 
used terms such as "zero modes" , " dispersionless modes" , 
and "partly flat bands" can be quantified. As shown, 
there is only one electron/hole state with zero energy 
at k y = it / a. It is characteristic for infinitely long 
graphene sheets with zigzag edges and belongs to the 
lowest /highest electron/hole band of in total 27V elec- 
tron/hole dispersion branches. The spectrum of any fi- 
nite stripe of zigzag GR contains no zero-energy levels at 
all. For example, in a 10 nm wide and 100 nm long zigzag 
graphene nanoribbon, the lowest/highest electron/hole 
level is E~ . , ast 10~ 97 (in this estimate, a = 1.42\/3 

rain/max ^ 1 v 

A 13]). The term "dispersionless mode or band" refers, 
in fact, to an exponential dispersion within an interval 
t(2N)- x w t/[4.7L(nm) + 2/3], up and down the Fermi 



level. The exponential dependence of edge-state energy 
on the wave vector changes from ~ exp(— ^/3Nq) to 
~ exp(2A ln(7r/3 — q) within an interval specified above. 
The precise quantification of the "flat part" of the low- 
est/highest conduction/ valence band is given by Eqs. © 
and M- 



2D GRAPHITE OR MACROMOLECULE 
MODEL? 

The above presented results refer to two models of 
graphene which are totally equivalent, as an object of 
Physics [lij], but require different formal treatments. The 
2D graphite model of graphene (a) (see 11, H, [3]) has 
been exploited by BF. This is a periodic structure of two 
nonequivalent atoms A and B. The alternative model 
(b) is equivalent to a 2D crystal with four nonequiva- 
lent atoms. Therefore, the corresponding disp ersion re- 
lations (or band structures) are different [l5(. As long 
as the bulk properties of infinite graphene are in focus, 
both models and both descriptions are equivalent. But 
once boundaries and/or size effects come into play, the 
graphene properties are more easily and naturally de- 
scribed starting from the macromolecule model rather 
than 2D graphite model. This is demonstrated by the 
present consideration and by obtaining accurate analyti- 
cal expressions of the band structure for achiral graphene 
ribbons and carbon nanotubes 



171 ] . The use of the Dirac 



equations demands some boundary conditions which can 
be formulated on the grounds of general requirements 
(such as current conservation, etc.) [3, [l||. However, 
from our point of view, there is a controversy in any 
attempt to rationalize boundary effects with the use of 
imposed boundary conditions. The use of the macro- 
molecule model is free of this and other difficulties of 
approximate descriptions. 

As already mentioned the possibility to use the Dirac 
equation for a description of graphene electronic proper- 
ties comes out from 2D graphite model of graphene. This 
means that the dispersion relation is obtained with the 
use of periodic boundary conditions for A and B sublat- 
tices. A commonly spread belief is that the dispersion 
relation is independent of the boundary conditions. But 
the case of graphene presents an exception. Modeled as 
a real macromolecule, the graphene dispersion relation 
has only two zero-point energies (instead of six at the 
corners of the hexagonal Brillouin zone) neither of which 
has k coordinates coinciding with the hexagon corners 
(To] ]. Therefore, references to K points in discussions of 
graphene electronics should be used with precaution. 

In summary, it is shown that for the larger part of the 
band of edge states in zigzag graphene ribbons, the de- 
scription given on the basis of the k-p approximation and 
corresponding Dirac equations agrees poorly with the ba- 
sic tight-binding model. In particular, such description 



4 



does not identify ak y — ±n as the only points, where 
the conduction and valence bands touch each other, and 
it does not reproduce correctly the exponential decrease 
(increase) of the edge-state electron (hole) energy to zero 
value. We have suggested an exact analytical description, 
based on the same model assumptions. For the larger 
part of the band in focus this alternative description is 
expressed in elementary functions. As a continuation of 
this work, accurate analytical expressions for low- (high-) 
lying electron (hole) bands have been obtained for zigzag 
and armchair graphene ribbons and carbon nanotubes 
171 ] . The use of the k-p approximation for these struc- 
tures becomes thus unnecessary. 
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